Mobility gap in fractional quantum Hall liquids: Effects of disorder and layer thickness 
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We study the behavior of two-dimensional electron gas in the fractional quantum Hall regime in 
the presence of finite layer thickness and correlated disordered potential. Generalizing the Chern 
number calculation to many-body systems, we determine the mobility gaps of fractional quantum 
Hall states based on the distribution of Chern numbers in a microscopic model. We find excellent 
agreement between experimentally measured activation gaps and our calculated mobility gaps, when 
combining the effects of both disordered potential and layer thickness. We clarify the difference 
between mobility gap and spectral gap of fractional quantum Hall states and explain the disorder- 
driven collapse of the gap and the subsequent transitions from the fractional quantum Hall states 
to insulator. 



I. INTRODUCTION 

One of the most remarkable properties of two- 
dimensional electron gas (2DEG) is the amazing preci- 
sion of the Hall resistivity quantization p xy = h/(ie 2 ) in 
a perpendicular high magnetic field at low temperatures, 
regardless of materials, geometries, impurities, and car- 
rier concentrations of experimental systems. This phe- 
nomenon is known as the integer quantum Hall effect 
(IQHE) for integer i, or as the fractional quantum Hall 
effect 2 (FQHE) for certain fractional values of i. At the 
Hall resistivity plateaus, the longitudinal resistivity p xx 
vanishes at zero temperatures, but has an Arrehenius- 
type temperature T dependence 

Pxx cx exp(-A/2fc B r) (1) 

where ks is the Boltzmann's constant liii 5 . The thermally 
activated behavior suggests that there is an activation 
gap A/2 in the excitation spectrum of each quantum Hall 
state. The gap has an origin of Landau level spacing 
in the IQHE 6 and a more profound origin of electron- 
electron interaction in the FQHEji in fact, the existence 
of the activation gap leads directly to the Hall resistivity 
quantization and the consequent incompressibility of the 
corresponding Hall liquid at low temperature. 

Theoretically, A in the FQHE is expected to be the 
creation energy of a pair of free quasielectron and quasi- 
hole. In a pure system, this is the asymptotic value of 
the excitation spectrum in the large momentum limit, 
found to be as large as OAe 2 /elsr^ where e is the dielec- 
tric constant and Ib = ific/ eB) 1 / 2 the magnetic length. 
However, experiments 3 ^^ found much smaller excitation 
gaps, presumably due to the reduction caused by the 
presence of disorder, the thickness of the 2DEG layer, 
and the mixing of Landau levels. Yoshioka™ combined 
the effects of layer thicknessii and Landau level mixing^ 



and obtained reasonable agreement with experimental re- 
sults in high-mobility systems at large enough magnetic 
fields £ 

However, without taking into account the effects of 
disorder, theoretical considerations cannot explain the 
vanishing activation energy below finite magnetic field 
(about 5 T)4£ Qualitatively, disorder broadens the 
quasielectron-quasihole excitation band, leading to a re- 
duction of the energy gap between the ground state 
and excited states^ MacDonald et nlML and Goldi^i 5 - 
considered the effects of disorder. Both theoretical ap- 
proaches contain adjustable parameters, and more im- 
portantly, fail to answer the following crucial questions: 
Which quasielectron-quasihole excitations are contribut- 
ing to the activated longitudinal resistivity? What is the 
nature of the activation gap? 

To answer these questions, let us look at the single- 
particle picture of the IQHE first. In a clean system with 
exactly n (labeled from to n—1) Landau levels filled, the 
activation gap is obviously the Landau level spacing hu> c , 
where lo c = eB/m*c is the cyclotron frequency. This 
involves the excitation of an electron in the (n — l)-th 
Landau level to the n-th Landau level, or equivalently, 
the excitation of a pair of electron in the n-th Landau 
level and hole in the (n — l)-th Landau level. In the pre- 
sense of disorder, each Landau level is broadened into a 
Landau band (with bandwidth 2r, say). Due to Ander- 
son localization, localized states exist in the tails of each 
band, and delocalized states only exist in the center of 
the band (where the mobility edge is). The excitation 
of an electron from a localized level to another does not 
contribute to the longitudinal resistivity; only excitations 
involving delocalized levels matter. Therefore, although 
the spectral gap in this case is reduced to (hcu c — 2T), the 
energy to excite a pair of free electron and free hole is still 
TiLj c - related to the mobility gap rather than the spec- 
tral gap of the system. In the FQHE, on the other hand, 
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we have quasielectron and quasihole excitations. Similar 
to electrons in the integer case, these quasiparticles can 
be trapped in their potential valleys and become local- 
ized, thus do not contribute to the longtidinal resistivity. 
Therefore, we need to find, for the FQHE, the mobility 
edge - the energy beyond which quasiparticle excitations 
are delocalized. 

In the noninteracting IQHE, the calculation of topo- 
logically invariant Chern numbers has been estab- 
lishe d 1 ^'^^?! 1 ?! 21 . 1 ! 2 ^ 22 ! 2 ? as a reliable way to obtain the 
Hall conductance, to measure the localization length crit- 
ical exponent, and to determine whether a single-particle 
state is localized or conducting - thus where the mobility 
edge is. Physically, the Chern number of a state is the 
(dimensionless) Hall conductance, which can be derived 
from the Kubo formular, averaged over boundary condi- 
tions of a finite system on a torus ^ In addition, it has 
an elegant geometric interpretation as the integral of the 
curvature of the quantum state in the parameter space 
spanned by two angular parameters (twisted boundary 
conditions) - the first Chern class of a U(l) principal fiber 
bundle on the torusi^25i 2 & Chern numbers are topologi- 
cally invariant under small perturbations of Hamiltonian, 
such as weak disorder, which allows the mobility gap to 
open for a finite range of magnetic field, leading to the 
plateau structure of the IQHE as long as the Fermi level 
lies in the mobility gap. 

However, the geometric interpretation of the Hall con- 
ductance leads to an apparent controversyi£*2£ in the 
FQHE. On the one hand, if a system exhibits the FQHE, 
the many-body ground state of the system must be 
degenerate on a toroidal geometry; otherwise, gauge- 
invariance arguments produce only integral Hall conduc- 
tance^ For a pure system at filling fraction v = p/q (p 
and q relatively prime to each other), the ground-state 
manifold is g-fold degenerate. Generically, these q de- 
generate states share a total Chern number p, regard- 
less of the Chern number carried by each state. Thus, 
the Hall conductance of the pure system is the average 
Chern number of these q states, fractionally quantized 
at ciff = pe 2 1 qh. However, numerical calculations show 
that impurities lift the degeneracy in a finite system^ 
implying that the Hall conductance would only be quan- 
tized to fractional values under superfluous conditions. 30 

Wen and Niu^i proposed that states of a system ex- 
hibiting the FQHE are topologically degenerate on a 
torus (in general, on high-genus Riemann surfaces) in the 
thermodynamic limit. The ground-state degeneracy is, 
in fact, a signature of the topological order of the bulk, 
invariant against weak but otherwise arbitrary pertur- 
bations (including symmetry-breaking impurity poten- 
tial)^ In a finite system, the quasidegeneracy of the 
states replaces the exact degeneracy, which can be re- 
covered in the thermodynamic limit. The topological 
degeneracy thus guarantees the Hall conductance to be 
quantized at fractional values in the fractional regime 
even in dirty systems, in which the topological Chern 
number for each many-body state is well-defined. Based 



on these ideas, we performed a numerical study of topo- 
logical Chern numbers for the v — \jZ FQHE in finite 
systems^ 3 " The mobility gap for the FQHE can be deter- 
mined from the distribution of the Chern numbers of the 
quasidegenerate many-body states. The results quanti- 
tatively explain the absence of the activation gap of the 
FQHE due to disorder at small magnetic fields, as well 
as the disorder-driven collapse of the gap and the subse- 
quent transition from the FQHE to insulator at higher 
fields, observed by various experiments^*^ 

In this paper, we further apply the method of Chern 
number calculation to study the effects of layer thick- 
ness of 2DEG and correlated potential on mobility gaps 
of FQHE systems. In Sec. |H] we introduce our micro- 
scopic model for 2DEG with disordered potential and 
finite layer thickness, and explain the Chern number cal- 
culation for FQHE systems. In Sec. IIIII we discuss the 
effects of layer thickness on the properties of 2DEG, in 
particular on mobility gaps, which we compare to the ac- 
tivation gaps measured by experiments. We then discuss 
the effects of correlated impurity potential on mobility 
gaps in Sec. lIVI before we summarize our results in Sec.lVl 

II. MODEL AND METHOD 

We consider a two-dimensional (2-D) polarized inter- 
acting electron system on an Ii x I2 rectangular area 
with generalized periodic boundary conditions (PBCs) 

r(Iy)*(r) = e*'$(r), (2) 

where T(Lj) is the magnetic translation operator and 
j = 1,2, representing x and y directions, respectively. In 
the presence of a strong magnetic field, one can project 
the Hamiltonian of the system onto the partially-filled, 
lowest Landau level. Therefore, we consider the following 
projected Hamiltonian in the presence of both Coulomb 
interaction and disorder 33 

i<j q 

+ £E eV/4c V tq - Rl > (3) 

i q 

where Ri is the guiding center coordinate of the i-th elec- 
tron, £/q is the impurity potential with wave vector q, 
and 

, , 27re 2 , . 
V(q) = —F(q) (4) 

is the Fourier transform of the electron-elctron interac- 
tion. The factor F(q) generalizes the Coulomb interac- 
tion to the case with finite electron layer thickness (to be 
explained in the next paragraph). We use the Gaussian 
white-noise potential generated according to the follow- 
ing correlation relation in q-space 

W 2 

(UM = (5) 
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which corresponds, in real space, to 
(U(r)U(r')) = W 2 5(r 



(6) 



where W is the strength of the disorder (in units of e 2 /e) 
and A — 2tt N s l 2 B is the area of the system. To study 
the effects of correlated potential, we also generate the 
Gaussian correlated random potential U(r) according to 



(UM = 

which leads to 

(U(r)U(r')) 
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(7) 



(8) 



where £ is the characteristic correlation length. Note that 
in the limit of £ — > 0, we recover the Gaussian white-noise 
potential. 

To describe the thickness of a quasi-two-dimensional 
electron system, we employ the Fang-Howard variational 
wave function 



(z) = (6 3 /2) 1/2 ^ b2/2 , 



(9) 



where b depends on the material properties and the car- 
rier density of the system. The parameter 6 _1 has the 
physical meaning of electron layer thickness, thus we will 
use, unless otherwise specified, the dimensionless layer 
thickness (3 = (bls)^ 1 hereafter. The experimentally in- 
teresting range of parameter is /? ~ 1. For the Fang- 
Howard wave function, the reduced Coulomb interaction 
in two dimensions becomes^ 



V(r) = (e 2 /e) / dqF(q)J (qr), 
Jo 



where 



F(q) = 1 



9_q W\ , «x 
8b 8b' 2 J \ b) 



(10) 



(11) 



and Jo is the Bessel function of zeroth order. This factor 
F(q), the same as in Eq. softens the bare Coulomb in- 
teraction between electrons, especially at short distances. 

We diagonalize the Hamiltonian with Lanczos algo- 
rithm and compute the Hall conductance an through the 
Chern number calculation, which offers an unambiguous 
criterion to distinguish between insulating and current 
carrying states in an interacting system42^ With a uni- 
tary transformation 



exp 



N r 
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(12) 



we can write the boundary-condition averaged Hall 
conductance for the fc-th many-body eigenstate as 
an(k) — C(k)e 2 /h, where the Chern number C(k) for 
the state is 



C(k) = — (b d9- 
4?r J t 



1 de 



1 de 



(13) 



Here, the closed path integral is carried out along the 
boundary T of the boundary condition space (the mag- 
netic Brillouin zone) < #i, 82 < 2n. C(k) is exactly the 
Berry phase (in units of 2-7r) accumulated for the state 
when the boundary conditions evolve along T. We sepa- 
rate the magnetic Brillouin zone into at least 25 meshes 
depending on the system size and calculate the sum of 
the Berry phase from each mesh. For the mesh sizes we 
choose, we find converged integer Chern numbers. We 
emphasize that throughout the paper we use the rectan- 
gular geometry, which facilitates the calculation of Chern 
numbers. A heuristic but qualitative discussion on the 
ground-state splitting and localized quasiparticle excita- 
tions also exists for a spherical geomotryi^ 



III. EFFECTS OF LAYER THICKNESS 

In this section, we discuss the effects due to the finite 
thickness of the 2DEG for v = 1/3. We use the Fang- 
Howard variational wave function, introduced in Sec. [nj 
to describe the electronic wave function in the perpen- 
dicular direction. We consider the Gaussian white-noise 
potential, and compare the disorder effects to the ideal 
two-dimensional casei^ In the following subsections, we 
present results on density-of-state, energy and energy 
split of ground states, spectral gap, distribution of Chern 
numbers, and mobility gap. These results are qualita- 
tively similar for cases with and without finite layer thick- 
ness. 



A. Density-of-states 

We diagonalize the system to obtain up to 30 lowest 
states for each sample for N e = 3-8. Figure^ shows the 
evolution of the many-body density of the 15 states as the 
layer thickness changes in the upper panel for weak dis- 
order W = 0.05 and system size N e = 6 . For W = 0.05 
and P = 0.0, a spectral gap is visible in the density-of- 
states with the low energy peak consisting of three nearly 
degenerate states. As /3 increases, the energies of states 
becomes closer to zero while the spectral gap narrows, 
reflecting the softening of the Coulomb interaction. For 
(3 = 2.0, even a weak W = 0.05 is sufficient to destroy 
the spectral gap completely. For comparison, Fig. 
shows the evolution of the density of the 15 states as W 
increases for zero layer thickness /3 = 0.0 in the lower 
panel. Disorder broadens the density of states and the 
low energy peak merges with other states at W > 0.11 
for jV e = 6. 



B. Ground-State Energy 

The effects of the finite layer thickness on the ground- 
state energy of the fractional quantum Hall effect has 
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FIG. 1: Density of the lowest 15 many-body states for N e — 6 
at v = 1/3. (a) W = 0.05 with /3 = 0.0, 0.5, 1.0, 1.5, and 2.0. 
(b) (3 = 0.0 with W = 0.05, 0.11, 0.16, 0.22, and 0.27. 

been considered by MacDonald and Aersi^ 5 . In the ab- 
sense of disorder, they found that the the Laughlin state 
energy at v = 1/3 for (3 = 1.0 reduces to 0.576 of its 
value for (3 = 0.0. Chakraborty2& found a similar reduc- 
tion ratio of 0.579 using the hypernetted-chain method. 

In the presence of disorder, we find a similar reduction 
in magnitude of the many-body ground-state energy Eq 
in finite systems. Figure^ shows the reduction of the 
ground-state energy as a function of (3 for N e = 6 at 
v = 1/3. Here, we defined the ground-state energy as 
the average energy of the lowest three levels, which are 
topologically degenerate in the thermodynamic limit. To 
compare our results with earlier works 3 ^^, we first add, 
to the many-body ground-state energy, a single-particle 
contribution from interactions of an electron and its im- 
ages due to periodic boundary conditions. This energy 
can be related to the Coulomb energy of the classical 
square Wigner crystal with finite layer thickness^ 5 - We 
then extrapolate our data for N e = 3 — 7 electrons to the 
thermodynamic limit according to 

E (N e ) = E {N e -> oo) + a /N e . (14) 

The results for ground-state energy per electron arc 
shown in Figure |2t>- By extrapolating the results to the 
clean limit by a quadratic fit, we obtain a reduction ra- 
tio of 0.572 for the ground-state energy from j3 = 0.0 to 
f3 = 1.0. The value is in good agreement with known 
results in the clean casei^ 5 ^ In the presence of the im- 



FIG. 2: (a) Ground-state energy (excluding interaction be- 
tween electrons and their images) as a function of (3 for N e = 6 
at v = 1/3. (b) Ground-state energy per electron, extrap- 
olated to the thermodynamic limit, as a function of (3 for 
v = 1/3. 



purity potential, disorder shifts the ground energy state 
down, leading to a slightly larger ratio (e.g., 0.596 for 
W = 0.11). 

In the current and earlier works 35,36 -, the ground state 
energy is negative, due to the assumption of a uniform 
neutralizing background charge, which cancels out the 
singular contribution of the Coulomb interaction among 
electrons. Mathematically, we substract the singular 
q = contribution from the electron-electron interac- 
tion, as illustrated in Appendix 1X1 Therefore, a reduc- 
tion of Eq in magnitude is, in fact, an increase of the 
ground state energy, as layer thickness increases. This 
is mainly because the corresponding reduction of the 
attractive interaction between electrons and the back- 
ground charge negates the contribution from the soft- 
ening of the electron-electron repulsion. Meanwhile, E 
decreases with increasing W because electrons take ad- 
vantage of the negative potential region in the ground 
state as shown in Fig. [5] 



C. Energy split of the ground states of finite 
systems 

Wen and Niu^i pointed out that the ground states of 
the fractional quantum Hall state are degenerate on a 
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FIG. 3: Energy split of ground states as a function of the 
square root of electron number (linear system size) on a 
semilog scale for various /3 at W = 0.05 and v = 1/3. 

torus even in the presence of disorder. However, this is 
strictly valid only in the thermodynamic limit. For a 
finite system on a torus of length L, the ground states 
have an energy split, or a bandwidth Eb, of order 

E b ^e- L ^' A ^\ (15) 

where A is the quasiparticle-quasihole pair creation en- 
ergy and m* is the effective mass of the quasiparticlc. 
The energy split comes from the tunneling process that a 
virtually created pair of quasiparticle and quasihole pro- 
pogate in opposite directions and annihilate on the other 
side of the torus. 

In figure |31 we plot the bandwidth Eb of the lowest 
three states, which are well separated from all higher 
levels at v = 1/3 and W = 0.05, as a function of ^/Ne, 
proportional to the linear system size. For various /3, the 
decrease of Eb with V^Ve ~ L is consistent with the ex- 
ponential decrease in Eq. (|15f) expected by the theory^ 
as found in Ref. |^ for zero layer thickness. Though, we 
cannot completely rule out a power-law decrease of Ej, 
based on these finite-size values. We find Eb increases 
with the layer thickness, again consistent with Eq. Ijl5|l 
assuming that the quasiparticle-quasihole pair creation 
energy decreases with increasing /?, as the Coulomb re- 
pulsion between electrons reduces. From (3 = 0.0 to 1.0, 
the slope in the semilog plot of figure|31reduces by roughly 
30%. The reduction is, however, smaller than that of ei- 
ther the spectral gap or the mobility gap which we will 
discuss in the following sections. 

We also note that a weak periodic potential, instead 
of disorder, gives rise to a similar exponential depen- 
dence of the splitting on system sizej^i implying, again, a 
degenerate ground-state manifold in the thermodynamic 
limit. The nontrivial ground-state degeneracy, existing 
on Riemann surfaces with genus 1 (torus) or greater—, 
is a signature of topological order— possessed by these 
states. The basic physics of the topological degeneracy 
and the finite-size splitting of the ground-state manifold 
in FQHE is the same as in chiral spin liquid, which is also 
a prototype for studying topological order and quantum 



error-correcting code— 

D. Spectral Gap 

Since the ground-state manifold on a torus consists of 
three quasidegenerate levels for v = 1/3, we define the 
spectral gap as the energy difference between the third 
and the fourth lowest energy states 

E S = E(4)-E(3). (16) 

At large disorder, E s simply becomes the energy level 
spacing, when the lowest three levels no longer form a 
separate band from higher levels. 

We extrapolate E s to the thermodynamic limit by fit- 
ting E s to 

E s {N e ) = E s {N e -> oo) + a s /N e . (17) 

Figure QJi shows E s (N e — > oo) as a function of f3 for var- 
ious W at v = 1/3. For small W, E s decreases as (3 
increases, or as the electron layer becomes thicker. For 
large W, E s remains smaller than or close to 0.003, re- 
flecting the closure of the spectral gap of the quantum 
Hall liquid. This finite but small residual value comes 
from the energy level spacing of the insulating phase, 
which may disappear with a more sophisticate definition 
of the spectral gap. In fact, for fixed layer thickness, E s 
first decreases with increasing disorder, then increases 
weakly with disorder after reaching its minimum. This 
minimum signals the critical disorder at which the spec- 
tral gap closes in the thermodynamic limit. 

E. Groups of Chern numbers and their statistics 

Due to the topological three-fold degeneracy of the 
each energy level in the thermodynamic limit for v = 1 /3 
quantum Hall liquids on a torus, we define the 7V s -th 
group of states as the (3N g — 2)-th, (3N g — l)-th, and 
(3iV g )-th states, and calculate the Chern number of such 
a group C(N g ) as the sum of the Chern numbers of the 
three states within the group. While for small enough 
disorder, the states within each group are degenerate in 
the thermodynamic limit, this no longer holds for disor- 
der strength large enough to distroy the fractional quan- 
tum Hall phase and the accompanying topological order. 

In earliar workjS we found the following properties 
for the statistics of the group Chern numbers for two- 
dimensional electrons at v = 1/3 with zero layer thick- 
ness. 

1. For weak disorder, the Chern number of the low- 
est group is always unity, carried by the three low- 
est states. This together with the fact that the 
three states become degenerate in the thermody- 
namic limit, as well as the fact that there is a finite 
spectral gap separating the three states to the rest, 
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FIG. 4: Spectral gap, extrapolated to the thermodynamic 
limit, (a) as a function of /3 for various W, and (b) as a func- 
tion of W for various j3 at v = 1/3. 

is the manifestation of the v — 1/3 fractional quan- 
tum Hall state on a torus. In the thermodynamic 
limit, each degenerate ground state carries a Hall 
conductance of e 2 /3h. This also holds for the Chern 
numers of upper groups for small enough disorder. 

2. Large enough disorder destroys the quantization of 
the Chern number of each group in individual sam- 
ples from upper groups down to the lowest one. For 
certain disorder, the probability of the group Chern 
number being unity decreases sharply around one 
group, the energy of which we can define as the 
mobility edge. We will discuss the procedure in 
greater detail later in Sec. IIII Fl 

3. The Chern number calculation is robust for small 
systems with as few as 5 electrons. The fluctuation 
of the group Chern number has little size depen- 
dence, which, therefore, can be used as an robust 
indicator of the degree of derealization. 

Figure 2 in Ref. l33l summarizes these properties, which 
can be used to determine the mobility gap. We repeat 
the calculation in the presence of finite la yer thickness. 
Figure shows, in analog to Fig. 2 in Ref. |33, the prob- 
ability distribution P(C) for group Chern number C of 
the lowest 5 groups of states in systems of 5-7 electrons 
for W = 0.05 and for (3 = 1.0 and 2.0. We find that the 
above-mentioned properties remain intact qualitatively 



FIG. 5: Probability distribution P(C) of total Chern number 
C for the lowest five groups of states in systems of 5-7 elec- 
trons for W — 0.05 and for /3 = 1.0 and 2.0. Energy increases 
from left to right. The arrow in each panel marks the group 
of states located at the mobility edge. 



in the presence of finite layer thickness. Quantitatively, 
the increasing layer thickness shifts the mobility edge to- 
ward the ground state. For example, for (3 = 2.0, the 
probility P{C — 1) of the lowest group is close to unity, 
suggesting that fractional quantum Hall phase survives 
at W = 0.05. Meanwhile, P(C = 1) of the second lowest 
group drops sharply to about 0.6, and may drop further 
lower for larger systems. This indicates the location of 
the mobility edge. 

We compare, in Fig. El the distribution P(C) for larger 
disorder W — 0.16 and 0.11 for (3 = 1.0 and 2.0, respec- 
tively, again for the lowest 5 groups in 5- to 7-electron 
systems. For the strong disorder, even the quantization 
of the lowest group no longer holds, which becomes much 
smaller than 1, indicating an insulating ground state in 
the thermodynamic limit. 



F. Mobility gap 

As discussed above, the fluctuation of the total Chern 
numbers in each group of states is an indication of the 
degree of derealization, which we can used to deter- 
mine the mobility edge. Similar as in Ref. I3H , we define 
Pcxt = 1 — P(C = 1). The value of P ox t is the probability 
of the breakdown of the Hall-conductance quantization 
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FIG. 6: Probability distribution P(C) of total Chern number 
C for the lowest five groups of states in systems of 5-7 elec- 
trons for W — 0.16 and 0.11 for f3 = 1.0 and 2.0, respectively. 
Energy increases from left to right. 



and thus a measure of the derealization of the charge ex- 
citations. This is analogous to the non- interacting IQHE 
case, where particle excitations from localized states to 
delocalized states lead to fluctuation of the total Chern 
number, or Hall conductance. 

For small disorder, P oxt remains close to for groups 
of states beyond the ground-state manifold. This re- 
flects that the mobility gap (which separates localized 
states from delocalized states) is different from the spec- 
tral gap (which separates the ground-state manifold to 
higher-energy states), as excitations across spectral gap 
may not lead to the fluctuations of the Hall conductance 
or contribute to the longitudinal conductance. Although 
the two gaps likely disappear simultaneously when dis- 
order is large enough to destroy the quantum Hall state. 
Therefore, we expect that P e xt rises sharply (probably 
abruptly) at the mobility gap in the thermodynamic 
limit, a signature imprinted in finite systems as well. 
Here, we find the sharpest jump of P cx t between one 
group and its lower neighboring group, and define the 
energy of the higher group as the mobility edge. We 
measure the mobility gap E m (N e ) from the ground-state 
energy to the mobility edge for the system of N e elec- 
trons. Since the states in each group are degenerate in 
the thermodynamic limit, we use the average energy in 
each group to calculate E m (N e ) to reduce finite-size fluc- 
tuations. We then extropolate E m to the thermodynamic 
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FIG. 7: (a) Mobility gap E m as a function of W for various 
layer thickness /3. E m is extrapolated from systems of N e = 
4-8 electrons to the limit 1/N e — > 0. (b) Dependence of E m 
on inverse mobility l//x° for various (3. The dashed line is 
converted from a fit to experimental data (taken from Ref.0). 
Here, we use an empirical mobility-density relation as well as 
a mobility-disorder relation in the Born approximation fi° = 
eTi 3 1 '(m* 2 W 2 ). The data points on j/-axes (for W = and 
1/fj, = 0) are the gaps for pure systems, whose calculations 
are illustrated in detail in Appendix El 



limit from iV e = 4-8 electrons. We plot the resulting E m 
as a function of W for various (3 in Fig. [7^,. The plot 
clearly demonstrates that finite layer thickness reduces 
the mobility gap as expected. In order to show the overall 
trend in the disorder dependence of the gap, particularly 
for weak disorder for various 2-D layer thicknesses, we 
include the gaps for pure systems on y-axes (for W = 
and 1/fi = 0). The calculations of the pure gaps, which 
differ from the Chern number calculations in disordered 
systems, are illustrated in detail in Appendix [BJ Com- 
paring Fig.E^ with Fig.^J}, we find that the mobility gap 
and the spectral gap differ significantly for small disorder 
and small layer thickness, although they both decrease 
with increasing disorder, as well as with increasing layer 
thickness. The two gaps appear to disappear at roughly 
the same disorder strength. 

In experiments, an energy gap A in the excitation 
spectrum of the correlated many-body ground state can 
be extracted from the temperature dependence of the 
magnetoresistivity, p xx cx exp(— A/2fcsT), where fc^ is 
the Boltzmann's constant and A/2 is often referred as 
the activation energji&ii^. This activation energy is re- 
lated to the mobility gap we calculated, which sepa- 
rates the ground state from its delocalized quasiparti- 
cle excitations. Boebinger et alX systematically studied 
the activation energy for v — 1/3, 2/3, 4/3, and 5/3 
and its dependence on sample mobility (i (an indica- 
tion of disorder) in a series of GaAs-AL-Gai-^As sam- 
ples. For a class of high-mobility (at that time - with 
mobility up to 10 6 cm 2 /V s) samples, they found that 
A ~ 0.049e 2 /d — 6 K, consistent with a simple phe- 
nomenological model^ that assumes a disorder-broadened 
excitation energy level with half- width T — 6K. Willett et 
al. studied a then untra-high-mobility sample and com- 
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pared the activation energy to theoretical results incor- 
porating finite layer thickness^ and Landau level mix- 
ingpifl*i£ While the agreement between experimental re- 
sults and theoretical caculations (in the absence of dis- 
order) is satisfactory for magnetic field stronger than 10 
T, they diverge significantly at smaller magnetic field, 
presumably due to disorder. Simple theoriesA&i4 with 
ad hoc treatment of disorder fail to account for all the 
discrepencies. 

With our numerical calculations which treat disorder 
and layer thickness on an equal footing, we can now at- 
tempt to compare our results to experimental ones quan- 
titatively. In experiments, the mobility /i dependence of 
A can then be extracted from the known dependence of 
fi on the electron density n of these samples, since n de- 
termines the magnetic field B (or lg) at the 1/3 family of 
fillings. For semi-quantitative comparison, we use a typ- 
ical dependence, \i — ^{n/no) 1 ^ , where /io = 600,000 
cm 2 /V s and uq = 1.5 x 10 n cm~ 2 , as extracted from 
Fig. 1 of Ref. 0. For comparison, we assume that in our 
simple disordered model both the (zero field) mobility 
and the (high field) mobility gap are dominated by short- 
range scatterers (appropriate for these then high-mobility 
samples). In the Born approximation (as derived in Ap- 
pendix Q, we have 



ehy(m* 2 W 2 ) 



(18) 



Figure compares this empirical formula of A(/x ) with 
the mobility gap we obtained in our calculation for var- 
ious p. We find that the experimental data fall nicely 
into the range of 0.5 < j3 < 1. This is fully expected 
in typical experiments: the variational parameter b^ 1 for 
a typical GaAs-AlGaAs sample with an electron density 
of No — 10 11 cm -2 is close to the magnetic length l B r^ 
which scales with B^ 1 / 2 (e.g. l B = 57 A for B = 20 T). 
In particular, Willett et al. obtained b^ 1 = 39 ± 1 A for 
their sample, giving (3 = (Mb)- 1 « 0.68 at B = 20 T. We 
do not include, however, the effects of Landau level mix- 
ing, which leads only to a small reduction of the mobility 
gap for clean samples 



IV. EFFECTS OF CORRELATED POTENTIAL 

In this section, we briefly discuss the effects of impu- 
rity potential with finite correlation length for v = 1/3. 
The correlated potential is, in particular, relevant to 
ultra-high-mobility 2DEG, such as GaAs based systems 
in which impurities are introduced remotely above the 
2DEG. The mobility of these systems is generally be- 
lieved to be limited by remote impurity scattering, rather 
than by short-range interface defects. To consider the 
correlated potential effects, we introduce the Gaussian 
correlated random potential 
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FIG. 8: (a) Mobility gap E m as a function of W for Gaussian 
correlated potential with £ = 1.0, in units of Ib, compared 
with Gaussian white-noise potential (£ = 0) . (b) Dependence 
of E m on inverse mobility for £ = 1.0 and 0.0. 



as described in detail in Sec. [H] Here, we restrict our- 
selves, for simplicity, to bare Coulomb interaction with- 
out considering finite layer thickness. We are mostly in- 
terested in the case of small correlation length £ ~ 1, in 
units of Ib, where quantitative results can be reached. 

The results we find, in particular for £ = 1, are qual- 
itatively similar to those with Gaussian white-noise po- 
tential.- 3 Figure compares the mobility gap E m as a 
function of W for Gaussian white-noise potential (£ = 0) 
and Gaussian correlated potential with £ = 1. The trend 
that increasing disorder strength destroys the mobility 
gap is generically the same for impurity potential with 
or without correlation. However, E m survives larger W 
(by a factor of roughly 50%) for £ = 1 than for £ = 0. 
This is not surprising since impurity potential correlation 
enhances the electron mobility, in particular for low elec- 
tron densities. As shown in Appendix [HI the mobility of 
the 2DEG, in the Born approximation, is enhanced by a 
factor of 



JL 



io(k 2 F e)-L (k 2 F ev 



(20) 



where Iq(x) is the zeroth-order modified Bessel function 
of the first kind and L (x) the zeroth-order modified 
Struve function. Here, hp — V^jm — \J~2v/l B is the 
Fermi vector of the polarized 2DEG. For £ = 1, the fac- 
tor is /i//i° — 2.556. In the long-range limit, one obtains 



/' 



(21) 



(19) 



Figure |SJd shows the mobility gap E m as a function 
of l//i for impurity potentials with and without corre- 
lation, as we have done here and in an earlier paper^ 
Interestingly, E m for £ = 1 scales back and lies roughly 
on top of the data for £ = 0. This suggests that as long 
as £ is not too large, the effect of the range of potential 
can be lumped into that of sample mobility; this makes 
comparisons between samples of different type, as well 
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as between theory and experiment more meaningful, be- 
cause the mobility is the directly measured character of a 
sample, while the details of disorder potential in general 
vary from sample to sample. 

We do not have, however, quantitative conclusion for 
vary large £, which is believed to be responsible for those 
ultra-high-mobility samples. Our finite-size calculations 
prevent us from considering large £ close to or even larger 
than the system size, at which stage the system prop- 
erties saturate. We would also need to calculate Chern 
numbers for a lot more eigenstates to obtain the mobility 
gap for the smooth potential, which is beyond our cur- 
rent computing capabilities. Nonetheless, it is expected 
that for those ultra-high-mobility samples, disorder plays 
a less important role, and thus the mobility gap depends 
less on the mobility of a sample. 5 



V. CONCLUSIONS AND DISCUSSIONS 

Semiconductor heterojunctions with modulation dop- 
ing is the de facto fabrication technique for high-mobility 
quasi-two-dimensional samples commonly used for the 
study of the fractional quantum Hall effects. The range 
and strength of the impurity potential and the finite layer 
thickness of the systems can significantly affect the trans- 
port properties of the system, such as the activation gap 
of the fractional quantum Hall liquids. In this paper, we 
have studied these issues using a microscopic disordered 
model of the fractional quantum Hall liquids in a toroidal 
geometry. With the help of the Chern number calcula- 
tion, which is capable of directly probing the localization 
properties of many-body states, we have determined the 
mobility edge of the Hall liquids based on the fluctuations 
of the Hall conductivity and studied the dependence of 
mobility gap on disorder strength, layer thickness, and 
potential correlation. 

Finite 2-D layer thickness has significant effects on the 
properties of the 2DEG in experimental samples. The 
spread of the electron wave function in the perpendic- 
ular direction softens the Coulomb interaction between 
electrons, reflected both in the density-of-states and in 
the ground state energy of the system. Consequently, 
the spectral gap, which separates ground state manifold 
from excited states, decreases with increasing layer thick- 
ness, as well as with increasing disorder strength. The 
mobility gap, associated with the fluctuation of topolog- 
ical Chern numbers - signaling delocalized excitations, 
also decreases with increasing layer thickness and disor- 
der. However, the two gaps are different by definition 
and, indeed, distinguishable in numerical calculations, in 
particular for small layer thickness and weak disorder. 
Putting reasonable experimental parameters, we found 
our results of the mobility gap are in excellent agree- 
ment with the activation gap measured by experiments, 
suggesting that the mobility gap is responsible for the ac- 
tivated behavior in longitudinal resistivity. Disorder and 
2-D layer thickness are the two dominant factors affecting 



the value of the gap. 

We also investigate the effect of the correlated disor- 
dered potential on the mobility gap. For small correla- 
tion length, we found an enhancement in the mobility 
gap compared to the case of uncorrelated potential for 
the same disorder strength. Such enhancement is consis- 
tent with the enhancement in mobility of the 2DEG in 
the presence of the correlated potential. As long as the 
correlation length is small, the effect of the correlated po- 
tential can be attributed to that of the sample mobility, 
which demonstrates it meaningful to compare the gap as 
a function of mobility among samples of different type 
and between theory and experiment. 

In this paper, we model the electron wave func- 
tion in the perpendicular direction by the variational 
Fang-Howard wave function. This is appropriate for 
modulation-doped GaAs-Al^Gai-^As heterojunctions, 
in which the Fang-Howard function is a very good ap- 
proximation to the numerical self-consistent ground state 
wave function^ For realistic values of layer thickness, 
the qualitative nature of the ground states and the low- 
energy excitations of the quasi-two-dimensional systems 
remains unchanged, in particular in the incompressible 
phases. 

Meanwhile, Shayegan et al£Q- studied the fractional 
quantum Hall effects in thick parabolic quantum wells. 
Due to the selective doping of Al and screening, elec- 
trons experience a flat potential and, therefore, the elec- 
tron density are roughly uniform in the quantum well. 
This results in a significantly larger layer thickness, which 
increases with electron areal density. In such a system, 
Shayegan et aZ42> observed a dramatic decrease in the ac- 
tivation gap with increasing layer thickness and thus the 
collapse of the fractional quantum Hall effect. He et al. 41 
studied the effects of layer thickness using a phenomeno- 
logical model potential (neglecting the accompanying ef- 
fects of Landau level mixing) and obtained qualitatively 
consistent results to the experimental measurements. Al- 
though we have not repeated our calculations with a dif- 
ferent wave function more suitable for a parabolic quan- 
tum well, we expect a similar trend of decreasing mobil- 
ity gap with increasing layer thickness. We would like to 
point out that Fig.[7|clearly demonstrates that increasing 
layer thickness can trigger a transition from a v = 1/3 
fractional quantum Hall liquid to an insulator for fixed 
disorder strength. The transition becomes easier to occur 
in the presence of larger disorder. 
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APPENDIX A: GROUND STATE ENERGY OF A 
TWO-DIMENSIONAL SQUARE WIGNER 
LATTICE WITH FINITE LAYER THICKNESS 

We calculate the ground-state energy of finite-size sys- 
tems with periodic boundary conditions. To compare 
our results with earlier works^^i, we must add, to our 
numerical ground-state energy, a single-electron contri- 
bution from the interaction of an electron and its im- 
ages due to periodic boundary conditions (see also Ap- 
pendix [0) . In an ideal two-dimensional square system 
of linear size L — \/2ttN s I b , this is simply the Madelung 
energy of a square lattice with lattice constant L and, for 
each electron, is 42,43 



For a = 1/2, the incomplete gamma function is also re- 
lated to the complementary error function 



erfc(z) = — = 

V7T 



dt e 



(A7) 



for z > 0. 

With the help of an integral transform, one can replace 
the slowly converging sum of Coulomb energy into two 
rapidly converging sums: one for the short-range part 
and the other the long-range, in the same spirit as in the 
original Ewald method^ In the following, the parameter 
that separates the two sums is denoted as yo, chosen to be 
tt/L 2 in our calculation. For convenience, one introduces 
the Jacobi 9 function 



e(z,x)= £ 



g27vmz ^ — rrrn X 



(A8) 



3.9e 2 



2L 



(Al) 



calculated first in the context of the Coulomb energy of 
the two-dimensional classical Wigner crystal^ The sum- 
mation is performed over lattice sites 1 = l\B.\ + for 
the primitive lattice vectors ai and &2 except 1 = and 
the 4> n (x) are the Misra functions 



+n ~— xt 

act e 



(A2) 



The factor of 1/2 in Eq. I|A1[) comes from the double 
counting of the electron-electron interaction. 

In the presence of finite layer thickness, the Coulomb 
energy of the quasi-two-dimensional classical Wigner 
crystal is also known for generic lattices i£ The wave 
function 4>{z) in the perpendicular direction enters 
through the following function 



f{y,b) = / dz / dz> 2 (z)0V)e 
Jo Jo 



2 ( -,'\r-y\ z - z '\ 



(A3) 



and for the Fang- Howard wave function 

(b(z) = (b 3 /2) 1 / 2 ze- bz/2 , (A4) 
one can obtain f(y, b) = f(b 2 /Ay) : where 



^=4- 8-2 + 2 r U'* 



(A5) 



Here, the incomplete gamma function T(a, x) is given by 



T(a, x) 



dtt° 



X a (f>a-l{x). 



(A6) 



which converges fast for not too small X. Here, we only 
review the results for a square lattice with lattice con- 
stant L for simplicity. To be specific, they are the q — > 
limit of Eqs. (16)-(21) in Ref. for the square lattice. 
The Coulomb energy per electron for the square lattice 
can be written as 



N7 1 E PP = NZ X E> + N~ 1 E< 



where 

El 
N e 

and 



2^ 



2L 2 



dyy- 1/2 f(y,b) 



i/n 



e z o 



L 2 y 



(A9) 
- 1 

(A10) 



~2L T 

„2 



dyy-^ 2 f(y,b) 

OO 

dyy- 3 / 2 f(y,b) 



L 2 y 



+N- i E':r( q )\ 



dyy-^ 2 f(y,b) 



q=0, 



(All) 



where E^° m (q) is the Coulomb energy corresponding to 
a homogeneous distribution 



N^EHriq) 



, dyy- 3 l 2 f{y,b)e- q /iy 

o 




15 



(A12) 



9=0 



Here the singular 1/q term has its origin in the lack of 
charge neutrality considered here. In fact, only nonsin- 
gular terms survive in the total Coulomb energy once 
neutralizing background charge (e.g., located uniformly 
at z — -- —bd) is present. The second term, as well as addi- 
tional nonsingular terms introduced by the interactions 
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between electrons and the background charge, depends 
only on "external" parameters, such as layer thickness b 
and location of the background charge bd- We neglected 
these terms since they have no effects on the results of 
the finite-size scaling in the 1/N e — * limit. 



p 


0.0 


0.5 


1.0 


1.5 


2.0 


A 


0.1037 


0.08410 


0.06729 


0.05656 


0.04940 



TABLE I: Gaps in pure FQHE systems at v = 1/3 for various 
2-D layer thickness. 



APPENDIX B: CALCULATIONS OF THE GAPS 
IN PURE SYSTEMS WITH LAYER THICKNESS 

In order to show the overall trend in the disorder de- 
pendence of the gap, particularly for weak disorder for 
various 2-D layer thicknesses, we presented the pure gaps 
(l//x = 0, W = 0) in Fig.0for v = 1/3. In this appendix 
we give the main details. These types of calculations in 
PBC geometry have so far been largely avoided for fear 
of strong finite size effects. These result from the inter- 
action of the quasiparticles with their images. Indeed 
such effects remain substantial (as large as 30% in some 
cases) for even the largest sizes in exact diagonalization 
studies. In order to remove these we follow the practice 
of not neutralizing the quasiparticle (qp) and quasiholc 
(qh) excitations 8 . This gives a positive contribution of 



(Bl) 



where e* = e/3 is the charge of the quasiparticle excita- 
tion, A is the area of the system, and 

V(r) = J d 2 q^-exp{iq-r}, (B2) 

where q = n\G\-\-niGi is the wave vector appropriate for 
the PBC unit cell, G's are the corresponding reciprocal 
lattice vectors, and n's are integers. We next subtract 
the repulsive interaction energy of a single quasiparticle 
with its images: 



— y 



1 



\X(l)\ 



(B3) 



where X(£) = i\L\ + £2^2, t's are integers and L's are 
the direct lattice vectors. These two terms together add 
to 



5E = (e*/e) 2 \S\ = |5|/9, 



(B4) 



where S is the classical ground state energy per electron 
(Madelung energy) of a Wigner crystal of electrons^^. 
Note that we are treating the quasiparticles as point ob- 
jects. This approximation only introduces another finite 
size effect since the size of the unit cell is several magnetic 
lengths whereas the substantial density variation of the 
quasi-particle excitations occurs over a magnetic length. 

Starting from the usual expression of the gap in terms 
of the ground state energies: 



A = E(v = l/3 + qp)+E(u = l/3 + qh) 
-2E(v = 1/3). 



With the above subtractions we obtain: 

A* = E*{v= 1/3 + qp) + E*(u = 1/3 + qh) 

-2E*(v= 1/3) + SE qp + 5E qh , (B6) 

where E* is the finite part of the energy without the 
Madelung term (E* = E — NS) (note S is negative). 
This contribution taken together for the 3 ground state 
energies in A* is a finite size effect as are both SE energies 
so they will disappear in the thermodynamic limit (A^ = 
Aqo). Therefore we need not correct any of these energies 
for finite layer thickness. We then extrapolate A* to the 
thermodynamic limit from its l/N dependence. We used 
the hexagonal unit cell for this part as it has the highest 
degree of symmetry. Table 1 gives the numerical values 
of the Aqo that we have obtained from TV = 7 — 10 size 
systems. In the same table we give the corresponding 
results for the finite layer systems. We have shown all 
these on the y-axis of Fig. [7\ (for 1/fj, = and W = 0). 



APPENDIX C: MOBILITY OF 
TWO-DIMENSIONAL NON-INTERACTING 
ELECTRONS IN RANDOM POTENTIAL 

In this appendix, we review the results of the mobil- 
ity of a two-dimensional non-interacting electron system 
with a random potential in standard perturbation the- 
ory^ and apply them to our model system. 

Consider spinless electrons in a quenched random po- 
tential U(r) with Gaussian correlation 



(U(t)U(/)) = 



w 2 

2< 2 



|r-r'| 2 /2£ 2 =WQt-t'\), (CI) 



where £ is the characteristic correlation length. Note that 
in the limit of £ — > 0, we recover 



(U(r)U(v')} = W 2 8{r - r') 



(C2) 



for a Gaussian white-noise potential. The total scattering 
rate in the Born approximation 



1 



27175 



Att 2 p 



--VF(2fc F sin^)(l -cos0) (C3) 
2tt 2 

drr[J 2 (k F r)-J 2 (k F r)]W(r), 



(B5) 



where p = m* /{2ttTl 2 ) is the density of states of free 
electrons and kp = y/^nn — V2v/Ib is the Fermi vector. 
W(k) is the Fourier transformation of W(r) 



(C4) 



W(k)= d 2 rW(r)exp(-ik-r). 
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For the Gaussian white noise potential, W(k) = W 2 , 
thus 



1 

Ttr 



m*W 2 



(C5) 



The mobility of the system with the short-range potential 
is, therefore, 



^~tr 

m* 



3h 3 



W 2 ' 



(C6) 



With simple algebras, we can rewrite the mobility fi in 
terms of the cyclotron energy ftio c , the Coulomb energy 
e 2 /e?s, and the magnetic length lg as 



— —it 



° 2 /eh 



1 

W 2 '' 



(C7) 



where W = eW/e 2 is the dimensionless disorder strength. 
For a long-range potential, kpS, ^> 1, one fmds2£ 



1 



(C8) 



In particular, for the Gaussian correlated potential, 

W{k) = W 2 e~ k ^ 2/2 , (C9) 
we can integrate Eq. (|C3() and obtain 
1 m*W 2 



Tti 



-e- k ^ [h{k%e) - L (k 2 F e)] , (CIO) 



where Iq(x) is the zeroth-order modified Bessel function 
of the first kind and Lq(x) the zeroth-order modified 
Struve function^ The mobility of the system with the 
Gaussian correlated potential is, therefore, 



H = n 



Io(kU 2 ) ~ L (kU 2 Y 
In the long-range limit, we obtain 



^{k F tf 



(Cll) 



(C12) 
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